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Abstract 

Modefs of superdense star clusters at the center of galaxies are investigated 
to see wliether such objects can be stable and long-lived based on evaporation 
and collision time-scales and stability criteria. We find that physically rea- 
sonable models of massive clusters of stellar remnants can exist with masses 
> 10^ Mq , which could simulate black holes at the center of galaxies with 
large M/L ratios and gas motions of order > 10^ km s^^. It follows that the 
evidence is not conclusive for massive dark objects at the center of galaxies 
being black holes. 

Subject headings: black hole physics - galaxies: individual (NGC 4258) - galaxies: nuclei 
masers 
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I. INTRODUCTION 



There has been a recent show of exuberance about the supporting evidence for the 
existence of ~ 10^ — 10^'^ Mq supermassive black holes at the center of galaxies (Kormendy 
1993; Kormendy & Richstone 1995 and references therein, Kormendy et al., 1997). However, 
as emphasized by Kormendy and Richstone, large M/ L ratios and gas motions of order > 10^ 
km do not provide a unique signature for supermassive black holes. Could these central 
dark objects in galaxies be massive clusters of stellar remnants, brown dwarfs, or halo dark 
matter? The arguments used to support the evidence that the massive dark objects are 
supermassive black holes are of an indirect astrophysical nature. The event horizons of the 
black holes would reside ~ 10^ — 10^ Schwarzschild radii below the spatial resolution of the 
Hubble space telescope. 

Maoz (1995) has used the measurement of the rotation curve of maser emission sources at 
the center of NGC 4258 (Miyoshi et al., 1995) to support the conclusion that this dynamical 
system cannot be a central cluster, unless the cluster consists of extremely dense objects 
with mass < 0.03 Mq, e.g., low-mass black holes or elementary particles. 

We shall show that it is possible to construct models of clusters with central masses in 
excess of 10^ Mq, which have reasonably long evaporation and coUision time-scales for such 
galaxies and can be expected to be stable objects. Miyoshi et al., obtained data using the 
Very Long Baseline Array, confirming that for NGC 4258 (M fti 4 x 10^ Mq) the rotation 
curve is Keplerian to a high precision. If the rotation is circular, then the mass interior to 
0".005 = 0.18 pc is M ~ 4 X 10''' Mq. This gets closer than is normally possible to a potential 
black hole at the center. 

Supermassive stars with arbitrarily large redshifts were first studied by Zel'dovich and 
Poduretz (1965) and subsequently by Ipser & Thorne (1968), Bisnovatyi-Kogan & Zel'dovich 
(1969), Ipser (1969, 1970), Fackerell (1970), Bisnovatyi-Kogan & Thorne (1970) and by 
Weinberg (1972). These studies were motivated by the suggestion of Hoyle and Fowler 
(1967) and Zapolsky (1968) that the emission lines of QSOs might come from the centers 
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of supermassive star clusters, whose gravitational fields would produce all or most of the 
observed redshift. Since then evidence strongly supports that QSOs are at cosmological 
distances and that the redshifts are due to the expansion of the universe. However, there is 
now a renewed interest to study dark massive objects at the center of galaxies, because of 
the possibility that they could be either black holes or supermassive clusters. 

II. EQUILIBRIUM CONFIGURATIONS OF SUPERMASSIVE CLUSTERS 

Let us give an overview of the supermassive cluster equilibrium problem. We shall 
assume that the core cluster is mainly supported by the pressure of radiation rather than 
of matter and that it is in convective equilibrium and has a uniform chemical composition. 
For radiation the energy density is e ~ 3p and the cluster can be described by a Newtonian 
polytrope with 



where m and n denote the mass of a star in the cluster and the number density of stars, 
respectively, p and p denote the mass density and pressure, respectively, and 1/(7 — 1) is 
a constant of proportionality. The condition of uniform entropy per star then yields the 
polytrope equation of state 



where X is a constant of proportionality. The standard polytrope equilibrium equation is 
(Chandrasekhar 1939): 



p — mn = (7 — 1) ^p. 



(1) 



(2) 




(3) 



where the radial variable r is related to ^ by 




(4) 



and 
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p = Kpie^'^^-^\ (5b) 

Boundary conditions on the Lane-Emden function 9 are given by = 1 and 9'q = 0. 

The total energy for the core Newtonian cluster is E — T + V, where T and V de- 
note the thermal energy and gravitational potential energy, respectively, and for Newtonian 
polytropes the total energy is E ^ -{3j - 4)GM^/(57 - 6)R. 

Since we have assumed that the supermassive cluster is mainly supported by thermal 
radiation, we can choose 7 ~ 4/3 and its mass is given by (Chandrasekhar 1939): 

K \3/2 



M = 25.3620 
and its radius by 



/A \'^/^ 



i? = 6.89685 (^—j po'^, (7) 

where po is the central density. 

The structure of the supermassive cluster can be described by a Newtonian polytrope 
with 7 ~ 4/3 (for NGC4258 GM/R ~ 1.6 x 10"^ which is locally Newtonian), but to settle 
the question of stability we need General Relativity (GR), for a cluster with 7 ~ 4/3 is 
sensitively balanced between stability and instability, so small effects due to GR and matter 
pressure must be accounted for, although they play little role in the structure calculations 
(Chandrasekhar 1964). 

We characterize the spherically symmetric cluster by a perfect fluid with the metric 
(Weinberg 1972): 

ds^ = B{r)df - A{r)dr^ - r\de^ + sin^ 0d(l)^), (8) 

where 



1 



2GM{r) 



r 



-1 



— [M(r') +47rr^p(r 



-1 



(9a) 
(9b) 



Moreover, we have 



M(r) = 47r f dr'r'^pir'). 
Jo 



(10) 



Outside the cluster p(r) and p(r) vanish and the metric is descibed the Schwarzschild solution 

(11) 



Bir) = A-\r) - 1 - '^^Mi^ forr > R. 



The general relativistic equation is given by the Oppenheimer-Volkoff equation (Oppen 
heimer & Volkoff 1939): 



p [r) 



^ _^ P(^)j N _^ 47rr^p(r)j r 2G'M(r) 
p(r)J L M(r) J L r 



(12) 



The thermal and gravitational energies are given by 



R 

T ^A-K I drr 



^ _ 2GM(r)l-V2 



(13a) 
(13b) 



where e(r) denotes the internal energy density. 
Expanding in powers of GM{r)/r we get 



■L 



R 

T = An I drr^ 

rR 

V — —An / drr 
Jo 



GM(r) 



+ ... 



GM{r) ^ SG^M^r) ^ 



p{r) 



(14a) 
(14b) 



r 2r^ 

A cluster described by a perfect fluid may pass from stability to instability with a radial 
normal mode at a value of the central density po for which the energy E and the number of 
stars N is stationary (Chandrasekhar 1964): 

dE 



dpo 

dN 

dpo 



0, 



= 0. 



(15a) 
(15b) 



The energy E has the approximate form 

rR fR i-R 

Ec^Anl drr^e{r)+AnG drrM{r)e{r) - AnG drrM{r)p{r) 
Jo Jo Jo 

-6nG^ I drM'^(r)p(r). 
Jo 

5 



(16) 



The total pressure is p — Pr+Pm where Pr and Pm denote the radiation and matter pressure, 
respectively, and the internal energy density is 



e = 3pr 



1 + 



(17) 



3(r-i). 

where F is the specific heat ratio of the matter and /3 = Pm/Pr- The total pressure is 

p = Pr{l + so that to first order in small P, the ratio of energy density to pressure gives 

(3r - 4) 



(18) 



3(r-i)' 

We now integrate by parts and to first order in /3 and employing the approximate Newtonian 
equation 

GM(r)p(r) , , 

p\r) ^ (19) 

to evaluate p,p and M(r), we get (Weinberg 1972): 

(3r-4) GM^ G^M3 

2{T-lf R • i^^) 

The stability criterion 

(?r 9po dr 
yields the minimum radius for stability 



(21) 



_ 20.4(r - 1) GM 

For a core cluster with M 3.6 x 10^ Mq, /3 ~ 0.1 and F ~ 5/3 we get the minimum radius 

i?^in = 2.3x lO'^pc, (23) 

which is deep inside a typical supermassive core residing in a galaxy (NGC4258). The ratio 
of -Rmin to the Schwarzschild radius Rs — 2GM is Rmin/ Rs — 6.8//?. 

Although we have used a perfect fiuid to investigate the stability properties of a super- 
massive cluster, we consider that it is a good approximation and that it is safe to conclude 
that a kinetic particle description of the supermassive cluster will yield similar results. 
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III. TIMESCALES OF EVAPORATION AND COLLISIONS 



The time-scale of evaporation of a bound system of objects with a single mass can be 



determined to be t, 



evap 



300trei<« (Spitzer & Thuan 1972; Spitzer & Hart 1971; Binney & 



Tremaine 1987) with 



t 



relax 



0.14A^ 



\gm) ' 



(24) 



.ln(0.4A^). 

where N is the number of objects in the cluster, M is the cluster mass, and R1/2 is the 
cluster radius within which lies half of the cluster's mass. 

Let us assume that all the mass in the dense cluster is confined within a spherically 
symmetric core, mantle and halo. The core is described by (Bisnovatyi-Kogan & Thorne 
1970): 

27rGpor^- 



Pc{r) = Po 



37 



(25) 



where po is the central density. This core is joined smoothly to a mantle with the density 
profile 

/ 1 



Pm(r) = 



7 



(26) 



,l + 67 + 7W27rG'r2 
with the join point in a region just outside the core radius. 

Outside the mantle a Newtonian envelope is constructed which is assumed to be convectively 
stable and has a finite radius re at which p goes to zero in a polytropic fashion. 



Pe(r) oc (re-r)^ > 



(28) 



We now assume that almost all the mass is confined within the core and the mantle. 
Then, we have p ~ Pc + Pm and the mass is 



drr'^Pc{r) + / drr'^pm{r) 



Y^'cPo + 



Sn^Gpl ,27, 

^^rc + -^(rm-rc). (29) 



From (^) and with 7 < 1 and Vm >> Tc, it follows that we can ignore the first two terms 
on the right-hand side of (pQ]). The evaporation time for R1/2 ~ 0.2 pc, M ~ 3.6 x 10'' 



and a dense cluster consisting of neutron stars with m ~ 1.4 Mq is 



evap 



lOGyr, 



(30) 



which is an adequate lifetime for the superdense cluster. 

Maoz (1995) has imposed constraints on the mass distribution of a dense stellar cluster 
by using the observational findings for NGC 4258. The high- velocity maser emission data 
obtained for NGC 4258 describe a nearly planar structure, and the velocity decreases from 
fin = 1080 ± 2kms^^ at a distance r^^ = 0.13 pc to fout = 770 ± 2kms~^ at a distance 
'"out = 0.25 pc (Miyoshi et al., 1995). Assuming circular motion and a perfectly planar 
disk, it was found that it can be fitted very well by a Keplerian relation. The systematic 
deviation of the velocity profile from a Keplerian relation is Av < 3 km (Maoz 1995) or 
a fractional deviation of Av/v < 4 x 10^^, where v is the average rotational velocity. Maoz 
assumed that the entire mass is within a radius Tin with a mass density profile described by 
a Plummer model (Binney & Tremaine 1987): 



The ratio of the cluster mass enclosed between the spheres of radii r^^ and rout to its mass 
within Tin is given by 



Then, for 5kcp ~ 0.01 and solving for Tc using (^Tj) and (p2D, it follows that < 0.012 
pc and with R1/2 = l-3rc and ^ Tc, we get for a hypothetical cluster of neutron stars 
(m ~ 1.4 Mq) at the center of NGC 4258 the evaporation time tevap ~ 10^ yr. Since this is a 
period of time much shorter than the age of the galaxy, Maoz ruled out the possibility that 
NGC 4258 is a cluster of stars with mass ~ 1.4 Mq. A cluster of stars with mass ~ 0.03 Mq 
would yield tevap ~ 6 Gyr, which would not be ruled out, but would be difficult to reconcile 




(31) 



"Kep 

= [M(< Tout) - M{< rin)]/M(< Tin). 



(32) 
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with collision timescales, unless the objects are extremely dense, e.g., light black holes or 
elementary particles, which are difficult to reconcile with any known theory of structure 
formation or stellar evolution. 

Consider now our model of a superdense cluster of stars with a mass profile given by 
(PSp. Assuming that r^. <^ ~ we get 

n„ < (33) 

-1- + f Kep 

and for the observational value ^Kep ~ 0.01 we obtain rin ~ rout ~ ~ 0.2 pc, which 
yields an evaporation time given by ( PP] ) consistent with the age of the galaxy NGC 4258. 
The structure of the supermassive cluster at the centre of the galaxy consists of a massive 
core and mantle with a tenuous gas envelope. We can conclude that such a compact object 
should not produce a deviation of the Keplerian rotation curve which exceeds the observed 
value Af < 3kms~^. 

Let us now consider the physical collision time-scale. This can be estimated from the 
formula (Binney & Tremaine 1987): 



coll 



(34) 



where n is the number density of stars, is the radius of the star, m is the mass of the star, 
and a is the velocity dispersion. For zero-temperature brown dwarfs and low-mass stars the 
mass-radius relation can be taken to be (Zapolsky & Salpeter 1969; Stevenson 1991): 



2.2 X 10^ 



/ m 







-1/3 



1 + 



m 



0.0032Mf^ 



-1/2' 



-4/3 



cm. (35) 



From (1^) and (^, we can estimate the core density to be 



12.93M 

For R ~ 0.2 pc, M ~ 3.6 x lO^M0,m ~ 1.4M0,r, ~ 1.97 x 10^ cm and a ~ 1500kms-\ 
we get po ~ 5.9 x lO^^Mgpc"^, n = po/m ~ 4.2 x 10^°pc~^ and 

teoii ~ 1 Gyr, (37) 
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which is long enough to offset a rapid evolution of the massive stellar cluster through co- 
alescence of stars. Of course, as we decrease the radius R of the massive core, then the 
time-scale of physical collisions of stars will decrease and lead to an unstable configuration. 

In a recent article (Moffat 1997) the fate of a dense cluster of stars which is undergoing a 
final stage of gravithermal catastrophe (Lynden-Bell & Wood 1968) was analysed. Since this 
phase of the evolution of a superdense core of stars is far from thermodynamic equilibrium, 
nonlinear cooperative contributions are expected to be important in the transport equations 
describing the last stage of evolution. It was found that such nonlinear contributions can 
prevent the core redshift from increasing without limits as the core becomes increasingly 
dense, preventing the collapse to a black hole. In particular, the redshift can remain less 
than the critical value for relativistic collapse, resulting in a stable, massive dark object at 
the center of a galaxy with a Newtonian core, mantle and thin halo. 

IV. CONCLUSIONS 

We see that it is not possible to rule out the hypothesis that NGC 4258 or other su- 
permassive galaxies are dense clusters of star-like objects with a mass m ~ 1 M©, since 
the structure of the theoretical mass profile of such superdense clusters is not known with 
certainty. Therefore, we must conclude that there is presently no conclusive evidence for 
the existence of black holes at the center of galaxies such as M31, M32, M87, NGC 4594, 
NGC 4258 and other potential black hole candidates. However, if these dark massive ob- 
jects with masses in excess of 10^ Mq are superdense clusters of stars with relatively large 
central redshifts, then they would be of considerable theoretical interest to the astrophysics 
community. 

Since it is difficult in the forseeable future to obtain an observational spatial resolution 
less than 10"^ Schwarzschild radii, it is not clear how the dark object at the center of galaxies 
can be proved to be a black hole. Black holes may exist at the center of galaxies but it 
cannot be claimed without further conclusive evidence that they have been been detected 
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by current observational data. 
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